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Abstract 

We propose a new model of leptons which corrects the LO predictions of Golden 
Ratio mixing via a minimal NLO Majorana mass correction which completely 
breaks the original Klein symmetry of the neutrino mass matrix. The minimal 
nature of the NLO correction leads to a restricted and correlated range of the 
mixing angles allowing agreement within the one sigma range of recent global fits 
following the reactor angle measurement by Daya Bay and RENO. The minimal 
NLO correction also preserves the LO inverse neutrino mass sum rule leading to a 
neutrino mass spectrum that extends into the quasi-degenerate region allowing the 
model to be accessible to the current and future neutrinoless double beta decay 
experiments. 
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1 Introduction 



It is an interesting feature of neutrino physics that two of the lepton mixing angles, the 
atmospheric angle ^23 and the solar angle 612, are both rather large jlj. Until recently 
the remaining reactor angle, ^13, was unmeasured. Direct evidence for 613 was first 
provided by T2K, MINOS and Double Chooz [2}1]. Recently Daya Bay [5], RENO [B], 
and Double Chooz [7] Collaborations have measured sin^(2^i3): 

Daya Bay: sin2(2^i3) = 0.089 ± 0.011(stat.) ± 0.005(syst.) , 
RENO: sin2(2^i3) = 0.113 ± 0.013(stat.) ± 0.019(syst.) , (1) 

Double Chooz: sin2(2^i3) = 0.109 ± 0.030(stat.) ± 0.025(syst.) . 

This measurement excludes the tri-bimaximal (TB) lepton mixing pattern |8j in 
which the atmospheric angle is maximal, the reactor angle vanishes, and the solar mixing 
angle is approximately 35.3°. When comparing global fits to TB mixing it is convenient 
to express the solar, atmospheric and reactor angles in terms of deviation parameters 
(s, a and r) from TB mixing [H]: 

1 1 r 

sin 6^12 = -^(1 + s), sin6'23 = ^(1 + a), sin 6^13 = (2) 

From the global fits in [T0Hl2] for illustrative purposes one may consider the represen- 
tative Icr ranges for the TB deviation parameters: 

s= -0.03 ±0.03, a = -0.10 ±0.03, r = 0.22 ± 0.01, (3) 

assuming a normal neutrino mass ordering. As well as showing that TB is excluded 
by the reactor angle being non-zero, Eq. ([3]) indicates a preference for the atmospheric 
angle to be below its maximal value and also a slight preference for the solar angle to 
be below its tri-maximal value. It seems that TB mixing no longer holds the exalted 
position that it did before, and perhaps now is an opportune moment to consider other 
mixing patterns that have been proposed but which have so far been somewhat eclipsed 
by TB mixing. 

An interesting alternative to TB mixing is the Golden Ratio (GR) [T314T5] mixing 
pattern: 



75 ^0, (4) 



/ 

Ugr = 



so named because it involves (f)g = , the famed Golden Ratio which ancient Greeks 
thought aesthetically pleasing. Notice that we have reported the above lepton mixing 
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matrix in the Particle Data Group (PDG) convention for Upmns II] where Pq is the 
matrix of Majorana phases. 

The history of the GR's possible role in lepton mixing beg footnote [13]. 

Several years later, this idea was applied [H] in the context of a non-diagonal charged 
lepton basis by observing that 6*12 = tan~^(l/0g) pa 31.7° was a good leading order 
(LO) prediction for the solar neutrino mixing angle. In addition to this, Ref. [H] con- 
jectured the possible connection of this prediction to the group A^. With this in mind, 
the authors of Ref. [I5] sought to elucidate the group theory of A5 and generate the 
aforementioned prediction for the solar neutrino angle prediction in the context of a 
non-dynamical flavour modelQ Shortly after this, it was found [T7j that there was an- 
other possible prediction for the solar neutrino mixing angle involving the Golden Ratio, 
612 = cos~^(0g/2) = 36°. But instead of using ^5, the dihedral group Diq was utilised to 
dynamically generate this prediction. For two years, the idea of the Golden Ratio's possi- 
ble implication in neutrino mixing lay seemingly dormantH, until a dynamical model 
was constructed [19] in a basis in which the charged lepton mass matrix is diagonal. 
Shortly afterwards, it was shown [20] that the prediction of 9^ = tan~^(l/0c,) ^ 31.7° 
was minimally realised in A^, when assuming a diagonal charged lepton basis. Yet a 
group in which the prediction 6^12 = cos~^(0g/2) was realised under the same assump- 
tions, could not be found. The same work [20] also contained two dynamical models 
which predicted Golden Ratio mixing. It should be noted that as all of this work using 
Ar, to construct models which predicted Golden Ratio lepton mixing was being done, 
was also being used to construct a four family lepton model [21], and its double cover, 
A'^, to construct a four family model incorporating quarks [22] ■ In addition to this, A'^ 
was used to construct a flavour model explaining cosmic-ray anomalies [23]. 

Unfortunately the measurement of the reactor angle also excludes the GR lepton 
mixing pattern in Eq.(jl]) in which the atmospheric angle is maximal, the reactor angle 
vanishes, and the solar mixing angle is given by 6*12 = tan~^(l/0g) ^ 31.7°, correspond- 
ing to r = a = and s ~ —0.09. These all lie outside the la ranges in Eq.([n]), especially 
the reactor parameter r, but also the atmospheric deviation parameter a, and even the 
solar parameter s is too negative. In the face of this disagreement, we shall interpret the 
GR prediction as the LO prediction and then try to achieve consistency with the data at 
next-to-leading order (NLO). This is analogous to a recent strategy that was developed 
for dealing with the TB prediction in the light of a non-zero reactor angle |24j, where 
the LO Klein symmetry based on the 6*4 generators S, U was broken at NLO but the 
generator S was preserved in order to maintain the successful tri-maximal solar pre- 
diction s = 0, leading to an atmospheric sum rule relation a = — (r/2) cos 5. However, 
in the present case, since we also wish to correct the solar angle prediction, we do not 

^It is important to note that Ref. 16 laid a large portion of the group theoretics of ^5 used in the 
exploration of Ref. [T^ . 

^With the exception of Ref. 18 constructing a quark model by extending A5 to its double cover, 
A5, with their golden ratio prediction of Ref. [15] in mind. 
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wish to preserve any of the original LO Klein symmetry, and so we shall break both 
S and U at NLO. For a general NLO correction, this would imply arbitrary values for 
the deviation parameters r, s, a. However we shall consider a minimal NLO correction 
leading to restricted and correlated ranges of these parameters. 

In the present paper, then, we propose a new model of leptons which corrects 
the LO prediction of GR mixing via a minimal NLO Majorana mass correction that 
completely breaks the original Klein symmetry generators, S and U, thereby correct- 
ing all three mixing angles and allowing agreement to be achieved with the global fits. 
Although there is no remaining mixing angle sum rule prediction, the minimal nature 
of the assumed NLO correction leads to a restricted and correlated range of the result- 
ing deviation parameters r, s, a which encompasses the la ranges of these parameters 
from the recent global fits. The minimal NLO Majorana mass correction also has the 
feature that it preserves a LO inverse neutrino mass sum rule even at NLO, leading 
to a quasi-degenerate neutrino mass spectrum. This, in turn, severely constrains the 
parameter space of neutrinoless double beta decay, allowing the model to be tested by 
next generation neutrinoless double beta decay experiments. 

The layout of the remainder of the paper is as follows. In Section [21 the Golden 
Lepton Flavour Model is constructed by defining its fields and transformation proper- 
ties under as well as an additional U{1) symmetry. The resulting LO and NLO mass 
matrices are constructed after electroweak and flavour symmetry breaking and diago- 
nalised to reveal the lepton mass and mixing parameters. In Section [3l the vacuum 
alignment of the model is analysed by explicit construction and minimisation of the 
model's flavon potential. Section H] contains a detailed analysis of the phenomenological 
implications of the NLO correction on the TB deviations parameters (r, s, and a) and 
the effective Majorana mass scale of neutrinoless double beta decay. In Section [5l the 
discussion of the A^ x U{1) model is concluded. The relevant group theory of A^ can 
be found in Appendix [Xj and a basis for the contractions of the NLO correction to GR 
mixing in Appendix [Bl In Appendix [21 the breaking of the low-energy Golden Ratio 
Klein symmetry is discussed. 

2 The Model 

2.1 Fields, Symmetries, and Yukawa Superpotentials 

In this section, we present an A^ model of leptons as a framework in which to study 
the generation of a non-zero ^13 from A^. We begin this discussion by noting that the 
left-handed lepton doublets L = {le,lfj,,lr), the right-handed charged lepton singlets 
^ = lAii ^r) 1 the right-handed neutrinos = ((^'H)e, ('^K)r) all trans- 

form under the 3-dimensional irreducible representation of A^. Furthermore, the up- 
and down-type Higgs doublet fields hu and are assumed to be blind to the flavour 
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symmetry. As a result of this minimal field content and its transformation properties 
h^LN, NN, and h^LE are the only allowed terms in the Yukawa superpotential be- 
cause 3 (g) 3 = Is + 3a + 5s contains a singlet (see Appendix |A] for the group theory of 
A^). However, these terms lead to undesirable phenomenological predictions, namely 
degenerate charged lepton and neutrino masses as well as no leptonic mixing. 

To fix the problematic predictions of the simplest LO terms {huLN, NN, and hdLE), 
it is necessary to introduce additional scalar fields {i.e. flavon fields) which will couple to 
the existing matter fields, as well as an additional U{1) symmetry to forbid problematic 
operators which lead to un-phenomenological results. The U{1) symmetry will not be 
gauged, in order to avoid constraints associated with anomaly cancellation. Yet, it 
will be spontaneously broken by the flavon fields acquiring vacuum expectation values 
(VEVs). In general, this will lead to massless Goldstone bosons, unless the symmetry is 
also explicitly broken. Therefore, the U{1) symmetry is assumed to be explicitly broken 
in the hidden sector of the theory, so that that Goldstone bosons become Pseudo- 
Goldstone bosons with mass around 1 TeV. The additional fields, as well as the original 
matter fields, and their transformation properties under As x U{1) can be found in 
Table [TJ Using these transformation properties, the Yukawa operators invariant under 
the Ar, X f/(l) symmetry can be constructed. 



Field 


L 


E 


N 


hu,d 


X 
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A 




A5 


3 


3 


3 


1 


5 


5 


1 


3 


5 


m 


-y 


+ y 


y 





— 5x 


— 3x 


—X 


-y 


-21/ 



Table 1: The field content of the x U{1) model where x and y are integers. 



The construction of the A^ x f/(l) invariant Yukawa superpotential is begun in the 
neutrino sector. Using the fields and transformation properties defined in Table [H the 
neutrino Yukawa superpotential may be schematically written to NLO as 

= yLNK + y2NNip + AM/,, (5) 

where /SWy denotes the NLO correction to the superpotential which will generate a 
nonzero 6^13 and a correction to 6*23 and O12] it will be discussed more fully in Section [231 
Proceeding to the charged lepton sector, the superpotential resulting in the generation 
of the charged lepton masses and mixings can be schematically expressed as 

Wi = ^LExhd + ^^LEO^ha + j^LEO'h,. (6) 

With the charged lepton and neutrino superpotentials in hand, we now turn to the 
construction of the charged lepton and neutrino mass matrices after electroweak and 
flavour symmetry breaking. 
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Flavon VEV 


VEV alignment 


(0) 


(0,0,0,t;<^,0)^ 


(x) 


(0,0,^;^,0,0f 


(A) 


(0,0, aaf 




(a/I (^2 + fs) , -"2, -"3, ^^3, V2V 



Table 2: The vacuum alignments of the flavons used in the model in terms of complex VEVs. 
Justification of the flavon alignments follows from a detailed discussion of the minimisation of 
the flavon potential in Section [3l 



2.2 The Charged Lepton Mass Matrix and its Predictions 

Using the superpotential in Eq. ([6]), the charged lepton mass matrix, after electroweak 
and flavour symmetry breaking, can be constructed. This is done by utilising the 
Clebsch-Gordan Coefficients of Appendix |X] to appropriately contract the flavon fields, 
thereby extracting the singlets from the product representations. Then, the VEVs of (j) 
and X given in Table [2], as well as (9) = vg are applied to reveal the charged lepton mass 
matrix, Mg, after electroweak and fiavour symmetry breaking, to be 

f^vl \ 

Me = vA VQj^v^v^ ^vl . (7) 

By assuming all fiavon field VEVs are of the same order, i.e. 7^ ~ e, and absorbing 
factors of -\/6 into and y^, a charged lepton hierarchy of e"^ : : 1 is obtained. With 
the added assumption that e ~ 0.15, a phenomenologically viable charged lepton mass 
spectrum is obtained. Furthermore, the associated charged lepton mixing matrix is the 
identity to O(e^). Therefore, charged lepton mixing may be neglected with respect to 
the mixing originating from the neutrino sector, as will be seen. 



2.3 The Leading Order Neutrino Mass Matrix 

In the case of (c.f. a similar situation in A4 [24j) a non-zero value for ^13 can be 
generated by not restricting a model to particular group representations but instead 
"trusting" the symmetry by adding all relevant irreducible representations. Since we 
choose to couple all of our fiavons to A^A^, where A^ ~ 3 and 3 ® 3 = Ig © 3a © 5s, 
(naively) three fiavons should be added transforming under the corresponding Ig, 3a, 
and 5s representations. Yet the 3a contained in 3 © 3 is antisymmetric (as denoted 
by the subscript "a"), and hence vanishes at LO from the tensor product of the two 
identical A^ fields. One could include the Is from 3 © 3 as well, as in Ref. [20], but 
it will not add any additional structure to the remaining 5s in 3 © 3; indeed, the 
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singlet contributes to the same Yukawa matrix elements as the first component of 
Therefore a trusting LO minimal GR model will contain only a single flavon ~ 5 which 

assumes a VEV (up to basis transformations) {(p) = {\J\{v2 + v^) ,V2jV^,v^,V2)'^ , 
where V2 and are left unspecified by A^. This VEV alignment is preserved by the 
five dimensional representations of the S and U generators in Appendix [Al and hence 
respects the corresponding x Klein symmetry of the LO neutrino Majorana mass 
matrix, 



M 



LO 
R 



1/2 



v 



{V2 + f 3) 
-\/3V2 
-V^V2 



It is straightforward to see that 



-\/3V2 



-\ iiv2 + V3) 



Diag(MfO,M2^^,M3^^) 



LO 



V^V2 
3 (f2 + ^^3; 



rLO\ 



\ 



(8) 



(9) 



where 



y2(^2(6(/>g-2)+4l>3) 

V6 



V6 



It is clear the above masses obey the sum rule 



rLO ^ ^LO 



M. 



LO 



Mi''' 



V3 



(10) 



Having constructed and diagonalised Mjf', it is trivial to construct the LO Dirac mass 
matrix, as it is just the Ig in 3 ® 3 resulting from the yLNhu operator in the super- 
potential (c.f. Eq. (jS])). After the spontaneous breaking of electroweak and flavour 
symmetries, the LO Dirac mass matrix has the form 



rriD = yvu 




(12) 



Applying the Seesaw Mechanism [25] to m^) and M^'^ generates the light neutrino mass 
matrix, 



M. 



LO 



/ V2-2V3 

' 3V2y2 



llvl 



2V3V2 - 4f ? 



v 



V2 

'2y2 

V2 

'2y2 



V2 
2^2 

(t)2-2l>3)(3l>2+2t>3) 

2V2{v2+iv3)y2 
13t;|+8«3^'2+4<;| 

6\/2(?J2+4u3)j/2 



V2 \ 
2y2 » 



l3vi+8v3V2+'iv'i 
6V2{v2+4:V3)y2 
(tl2-2l)3)(3l'2+2l)3) 

2v^(t)2+4t)3)j/2 



(13) 



^Also, the singlet does not break the S 01 U symmetries and therefore wiU not generate non-zero 
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that is diagonalised by 



GR 



P' 



( 



\ 




(14) 



where P' = Diag(l, 1, —1) is a matrix of unphysical phases and Pq = Diag(l, e~ , e 2 ) 
is the matrix of Majorana phases required to bring Ugr to the PDG convention [1] of 
Eq. dl]). However, to diagonahse M^*^ we must use the mixing matrix given in Eq. 
f ll4p and not Eq. as absorbing the unphysical phases into the charged leptons will 
interfere with the diagonalisation of M^*^. With this in mind, 



GR 




























(15) 



with complex masses 



m 



LO 



m 



LO 



LO 



(u2+4v3)j/2 ■ 



(16) 



Notice that these complex masses obey the inverse neutrino mass sum rule [20] 

rrJf mlf^ m\^' 

For the remainder of this work, it is useful to re-express the complex neutrino masses 
in terms of new parameters /3, ^, and b such that 



nil 



-2+4e«5' "^2 4e»'55-(^+2)' '"^ 2(1+46^''^) ' 

To arrive at the above forms of the complex neutrino masses, let f 2 = 
Then, define ^ = |f3|/|f2| and 5 = 6'i — 62. Then, it is clear /3 



:i8) 



\V2 



|e*^^ and 



^^3 



idz 



2/2 f 2 



Notice 



that and 6 are real parameters, but (3 is complex. Furthermore, it is interesting to note 
that the phase associated with P and is irrelevant when comparing to the experimentally 
measured parameters, |Amjjp. 

The re-parametrisation of the complex neutrino masses in terms f3, ^, and 6 con- 
cludes the construction of the LO lepton model. Notice that this model yields a 
zero reactor angle, maximal atmospheric mixing, and a solar mixing angle given by 
612 = tan^^(l/0g) ^ 31.7° which are all in conflict recent global fits. Therefore, it is 
necessary to consider the NLO corrections to the existing minimal LO GR model re- 
sulting from the addition of operators which will break the S and U generators of the 
low-energy x Z2 , providing a correction to the problematic leading order predictions. 
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2.4 The Neutrino Sector at NLO 



In order to generate a non-zero reactor mixing angle in the context of ^5, a flavon field 
must be added which breaks the U generator of the x Z2 Klein symmetry. To do 
this, an additional flavon field A transforming as a 3 under is introduced. Recall 
that a flavon field transforming as a 3 necessarily vanishes when coupled linearly to A^A^ 
due to the anti-symmetry of the 3a contained in 3 ® 3. Thus, we allow for a quadratic 
coupling of A to A^A^. As a result, the superpotential receives a next-to-leading-order 
correction oiQ 

AW. = ^ (NN),^ (AA),^ + ^ {{NN),^ {XX) . (19) 

This corrects the LO Majorana mass matrix in Eq. to yield 

Mr = ^ ^M^^ (20) 

where is defined in Eq. ([H]). Before proceeding to calculate AMfi, it is useful to note 
that the coefficient i/a will not enter the eigenvectors of the corrected/perturbed mass 
matrix. Therefore, it will not affect the corrected mixingH Thus, we choose an alignment 
for (A) = (01,02,03)^ such that the operator associated with i/a vanishes. This will 
require ai = and either 02 or 03 to be zero. We choose 02 = implying (A) = (0, 0, 03)"^. 
This choice will be justified in the detailed discussion of the minimisation of the flavon 
potential in Section [31 

Having assumed A develops VEV (A) = (0, 0, a^)'^, the correction AMr takes the 
rather simple form 

/ \ 

AMr = 6ycv^ , (21) 
\ / 

where ^4 = af/A. As discussed in Appendix [0 this correction turns out to be sufficient 
to break both the S and U Klein symmetry generators associated with the unsuccessful 
LO Golden Ratio predictions. 

Recall that the LO Golden Ratio mixing matrix is given by the usual GR form as in 
Eq. (jlj), up to charged lepton corrections (which are small as can be seen in Eq. ([7])), 
where Ugr has been brought into the PDG convention pQ. Notice that the P' of Eq. 
f|T^ has been absorbed into the charged lepton fields. In the presence of A, the PMNS 
mixing matrix is given as 

U'pMNs = Ugr + AU, (22) 

"^The pair of contractions given here form a basis for all possible contractions, see Appendix |B] 
Appendix [B] also makes clear the choice of notation for the Yukawa couplings jja and yc ■ 
^However, its existence will shift the neutrino masses depending on the alignment of (A). 
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where the prime on U'pj^^g indicates that it is not yet in standard PDG form due to 
the correction. Assuming 



(23) 



enables a perturbative calculation of the correction /\U . To first order in f4/f2,3, the 
corrections to the columns of the PMNS matrix are 



Af/ii 



e 

5' 



,75' 



Af/i 



12 



305 



2v^' 
Af/i3 ^ 18eV3L>^\ 

3 




(120, + 1 - 2ie'')D^\ 



Af/22 ~ 
Af/23 ~ 3ey6(l - 2ie'^)Dl\ 



At/31 ^ ev^(We'' - '/'s - 18)^2"', 
3 

Af/33 ~ 3eV6(l - 2ie'^)Dl\ 



Af/. 



32 




1 - 18(/., - 2ee*'^)D3-\ 



where 



Di = 272(4^^6^^^ + 2^e'^ - 11) 
D2 = 10 X 5'^\(f)g + 3 - 

= 10 X 5^/^(30g - 1 + 2^6*"^), and 

ycVi 
e = . 

V2y2 



(24) 

(25) 

(26) 
(27) 

(28) 

(29) 
(30) 

(31) 

(32) 

(33) 
(34) 
(35) 
(36) 



Therefore the correction to the PMNS matrix is 

AU = [AUa Af/,2 Af/,3) . (37) 

As noted in Section [H a general lepton mixing matrix can be expressed in terms of 
the TB deviation parameters r, s, and a. Doing so allows one to write the lepton mixing 
matrix as [9] 

1 „^ 



^^PMNS 



-;ig(l + .-a + re^^) ^^{\ - \s - a - \t e^') ^ 



1 + a 



P . (3J 
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where the diagonal matrix P = Diag(l, e~, e~) contains the usual two Majorana 
phases in the PDG convention pQ. The result of Eq. (155]) can be compared to the 
corrected GR matrix to show how the breaking of the S and U generators corrects the 
leading order prediction of r, s, and a. However before we make this comparison, phase 
conventions must be matched. Again, we use the freedom to multiply unphysical phases 
from the left since they may be absorbed into charged lepton rotations. Multiplying 
phases from the right corresponds to a redefinition of the Majorana phases. Once this 
redefinition has been made, the NLO PMNS matrix is in the PDG convention. Then, 
we may identify the deviation parameters r, a and s resulting from the corrected Ugr 
mixing matrix. They are 

r = y2|A[/i3|, (39) 

a = V2Re ( AU23 exp ( ^] ] , (40) 

AUuexp(^)), (41) 




where a2 and 0:3 are the redefined Majorana phases in the PDG convention [T]. 

We pause here to comment on the lack of a concrete prediction even though we have 
only introduced one parameter. Notice that the LO neutrino mass matrices given in 
Eqs. dH]) and f|T2|) respect Form Dominance [27] because in the basis that the Majorana 
mass matrix is diagonal, the columns of the Dirac mass matrix are proportional to 
the columns of Upmns- As such the GR mixing matrix of Eq. f lT^ contains no free 
parameters. By adding the correction, AM^ (c.f. Eq. fl^T]) ) to Mr, "one" additional 
free parameter is introduced {ycv^). However because of the way that the parameter 
enters, i.e. only in the (22) entry of Mr, the Majorana mass matrix can no longer be 
diagonalised with a transformation independent of parameters of the model, and the 
resulting Upmns will no longer be parameter free. Hence, at NLO the mixing angles 
will depend on parameters in the model, as Form Dominance has been broken. 

A similar (slightly easier) perturbative exercise may be performed to calculate the 
NLO corrections to the light neutrino masses. Performing this calculation reveals the 
light, complex masses at NLO to be 

^NLO ^ ^LO ^ 

^ ^ 2<t>V30[l-3(j}-2^ei^) 

^NLO ^ ^LO _ 9/30 ^42) 

2V3O(2-30+2?e»'5) ^ ' 

^NLO ^ ^LO 9/3 ^ 

^3 ~ A/t3 2y6(l+4?e»«)2^' 



where the are as defined in Eq. f[T6l) . Notice that the LO sum rule still exists at 
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NLO, I.e. 

+ - (43) 

This surprising result can be explained by first observing that the NLO correction to 
the heavy, right-handed neutrino masses will be given by the corresponding diagonal 
element of the matrix Uq^AMjiUgr- A calculation of these corrections reveals 

for the corrected heavy neutrino masses. Clearly M^^^ + M^^^ = M^^^, preserving 
the LO sum rule for the heavy right handed neutrino masses of Eq. f|TT]l . Furthermore, 
because of the form of rriD, it can be shown that. 



rrii 



vyMr\ (45) 



Hence, Mf^^ + M^^^ = M^^° implies l/mf^° + l/m^^^ = l/m^^^. In summary, 
due to the form of the correction to AM/j (c.f. Eq. fl2T]) ). the resulting matrix which 
governs the corrections to the heavy neutrino masses, Uq^AMrUgri also shares the 
same sum rule between its diagonal elements, allowing for the preservation of the sum 
rule for the heavy neutrino masses at NLO. Then, because of the form of mD (c.f. 
Eq. ffT^ ). it is possible to deduce that the heavy neutrino masses, Mj are inversely 
proportional to m^, preserving the light neutrino mass sum rule to 0(f4/f2,3)- 



3 Vacuum Alignment 



Field 


0' 


u 


0" 




x" 




A" 


e 




A, 


5 


4 


5 


3 


4 


3' 


1 


3' 


5 


U{1) 


— 6x 


z 


6x 


llx 


8x 


2y-z 


2y 


5x + y 


-2z 


U{1)r 








2 


2 


2 


2 


2 


2 


2 



Table 3: Auxiliary flavon and driving fields of the x U{1) model where x, y, and z are 
integers. 

Flavour models of this type must have the alignment of the VEVs of the flavon 
fields justified by minimising a flavon potential. Therefore, the explicit VEVs quoted in 
Section [2] must be derived from a flavon potential. Herein lies the goal of the present 
section. 

To properly align the VEVs of the flavon fields of the x f/(l) Golden Model, a 
set of auxiliary flavon fields will be added as well as a set of "driving fields". Driving 
fields are similar to flavons in that they are gauge singlets and transform in a nontrivial 
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way under x U{1). However, their difference becomes apparent when an additional 
U{l)ji symmetry is introducedjj Under this symmetry, we define all chiral superfields 
containining Standard Model fermions to have a U{1)r charge of +1. Then, all chiral 
superfields containing Higgs and flavon fields have a U{1)r charge of 0. Driving fields 
will have a t/(l)_R charge of +2. Since the superpotential carries a +2 U{1)r under 
these conventions, the driving fields may only couple linearly to fiavon fields. With the 
further assumption that the driving fields develop positive soft supersymmetric breaking 
scalar mass squared parameters at the symmetry breaking scale, the driving fields do 
not develop VEVs. Therefore, to minimise the fiavon potential it is only necessary to 
enforce that the F-terms of the driving fields vanish identically. These so-called F-term 
conditions will give rise to the vacuum alignments. The driving fields, additional fiavon 
fields, and their transformation properties under x f/(l) x f/(l)i? can be found in 
Table El 

To begin the minimisation of the fiavon potential, the first step is the construction of 
the the fiavon superpotential from the fields' transformation properties given in Tables 
[U and [31 A straightforward calculation shows the "LO" fiavon superpotential may be 
constructed from contractions of the operators 

000°, 9x(p', (P'XP'AXX', , u;^¥>o , AAA°, and Ax^"- (46) 

Because some of these operators have multiple independent contractions that can result 
from them, it is necessary to explicitly write the fiavon superpotential in terms of the 
relevant contractions. Doing this yields 

iy, = M(0y)^^+^?i ((00)5^^0°)^ +^?2((00)5,,>°)^ +93{0{x<P'),^)^ 

+9. m).j),^+9. m,y),^+9e m,J\^+97 ((^^)5.^°L ^^^^ 

+ ^78 ((c^<^)3' <^o)i^ + +99 ((AA)i^ A°)^^ + (710 ((Ax)3' • 

As previously discussed, the superpotential can be used to construct the driving fields' 
F-terms, that when vanish, minimise the fiavon potential. The vanishing F-term of p° 

^J7(l)fl is broken to R-parity when supersymmetry breaking terms are included. 
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and x° that provide the ahgnment for the VEVs of 0, 0', and x s-^e 



dp\ 


= 94 




- 201X3 + 203X4 + 02X5 


) = o 




dWd 


= 9a 




+ V30'lX5' 






= 94 


V 0</;2a1 


V 'JSfiA2 V z.</^5X3 ^ 


- V203X5) 


= 


dx\ 


i/5 




)5 - V3X204 + 4^3X30 


3 V '-'/V4y'2 


-J V ^a5V^1 I 




+ 96 




h + V3X204 - V^X402 


+ V2x50i'; 


= 


dWd 
dxl 


i/5 




)4 - V3X203 - V3X302 




'± V '-'a5V^5 / 




+ 96 


V2Xi04 


+ V3X203 - ^3X302 - 


V^X40l) = 


= 


dWd 
dxl 


= 95 




)3 + 4a/3X202 + 3^2X3^ 


f)i-V3 (x4 


05 + X504)) 




+ 96 


V2Xl03 


- V2X301 + V3 (X504 


- X405)) = 





dWd 
dxl 


= 95 




)2 + 3^2X201 - (X3< 


h - 4X404 


+ X503)) 




+ 96 


^-v^xk; 


62 + ^2X201 + (X50 


3 - X305)) 


= 



(48) 



which admit the solutions 

(X) = (0, 0, v^, 0, 0) , (0) = (0, 0, 0, v^, 0) and (0') = (0, v^,, v^,, 0, 0) . (49) 

Furthermore, it has been verified that 17 other group elements acting on the preceding 
VEV alignments also minimise the potential, providing a complete set of solutions|l| We 

^These solutions can be obtained by the action of T, T^, T^, T^, S, U, ST, TS, UT, TU, SU, T^S, 
T^U, UT^, T^S, T^ST^, or T^ST^ on the VEV ahgnments in Eq. dH]). 
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turn now to the F-term conditions resulting from (fP. These are found to be 
d<tP, 



dct>l 



M(t)\ + 2^1 {4)1 + 0205 - 20304) + 2^2 (0? " 20205 + 0304) + ^3^Xl = 0, 
M0'5 + 2^1 (0105 + 760204) + 92 (-40105 + V602) + 53^X5 = 0, 
M0; + gi (-40104 + V^05) + 2^2 (0104 + V60203) + ^3^X4 = 0, 
^03 + 91 (-40103 + V^(t>i) + 2^2 (0103 + V^0405) + ^3^X3 = 0, 
M0'2 + 2^1 (0102 + V60305) + 92 (-40102 + ^60^) + ^3^X2 = 0. 



(50) 



Applying the previous results for the alignments of (%), (0), and (0'), collapses these 
equations to the rather simple constraints 



Mv^i + g-sVgv^ = and Aff + V6g2vj 







> - - (51) 

which together imply = '^^^'^e'^x- Notice that since v^^ and 7^ 0, then vg 0. 
With the justification of the alignment of (%), (0), and {9), we turn ourselves to aligning 
the flavons associated with the neutrino sector. 

Aligning the flavons associated with the neutrino sector is begun by minimising the 
F-terms associated with ■0°. We flnd the following conditions on the alignment of the 
auxiliary field cu: 

(^V3u)iU)4 — 2^/3uJ2U)3j = 0, 
(2V2UJI - 2V2uJiUJs^ = 0, 
(2V2a;ia;2 - 2V2ujfj = 0, 
(2V2UJ3UJ4 - 2V2ujI^ = 0, 

(2^20;^ - 2V2U2U}4^ = 0. 









= 97 








= 97 






d^3 


= 97 


dWd 






= 97 








= 97 



(52) 



which have solutions of the form 

/ 1 \ 

1 



1 

VI/ 



p 



f P \ 

p2 

p3 

\ / 



/1\ 

1 
1 

\1/ 



(53) 
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in which p — and m,n — 0, 1, 2, 3, 4. For the rest of this work, we take m — n — 
such that (u) = Vi^{l, 1, 1, l)"^. 

With the ahgnment for (u) calculated, we turn now to the F-terms associated with 
(f^ and find the following F-term conditions: 



9s 



98 



- US (^-2(fi3UJi + (P2UJ2 + V6(fiiUJ3 - 3(/75a;4^ = 0, (54) 

= QS (3(p2i^i - V6ipiUJ2 - ip5i^3 + '2(p4UJ4^ = 0, 

By assuming the previous solution for {cu), the above equations reduce to 



dWr 



(y6(pi + (p2- 2</?3 - 3</?5^ = 0, (55) 



Since v,^ ^ 0, the above equations may be simplified to 



ip2 + 2(/?3 - 2(/74 - (/75 = 0, 

a/6v9i + - 2933 - 3</75 = 0, (56) 
V&^pi — 3(p2 - 2(^4 + (p5^0. 

These equations yield, up to action by T4, the solution 

{(P) = (^2 + -^3) , -^2, V3, V3, . (57) 

The remaining set of F-terms to analyse belong to A° and They are 



9A0 

0^2 



gg {Xl + 2A2A3) , 
gio (ySXiXi + A3X2 + A2X5) 
gio (^-^/2A2X3 + A1X4 - v^AaXs) 
g-io (^-V2X2X2 + A1X3 - V2A3X4) 
16 



(58) 



Upon applying the previously found alignment for (x) {i-e. (x) c>c (0,0, 1,0,0)"^), these 
F-terms vanish when (A) is aligned as 



the Z2 X Z2 breaking alignment for (A) used in Section 12.41 Before concluding the 
discussion of the flavon potential, we pause to comment on the generality of the A^xU (1) 
model's vacuum solutions. 

In the preceding section's discussion of the alignment of the VEVs of 4>, and (f)', 
it was noted that there exist 17 additional group elements that, when acting on the 
reported VEV alignments, will also satisfy the relevant F-term conditions. As such, 
the 18 group elements represent the complete set of 18 solutions, up to proportional- 
ity. Without loss of generality, we can choose the solution given in Eq. (H9l) . Having 
made this choice, the VEVs of 6 and A are determined uniquely. Next, the VEV of u 
may be brought to the form (1, 1, 1, 1) by action of T transformations, which in turn 
uniquely determines the alignment of (ip). Yet, having done this we must check that T 
transformations do not change any of the previously aligned solutions. Action by the T 
generator on the VEVs 4>, and (p' introduces irrelevant factors of p on all of the VEVs 
with the exception of {</)'), where action by T introduces a relative phase between the 
two non-vanishing components. One could argue that because this is an auxiliary field, 
only introduced to aid in the minimisation of the potential, there is no problem since 
all five solution yield the same physical results for lepton masses and mixings. However 
before this claim is made, it is insightful to look at the previous argument in "reverse" . 

Because the alignment of (u) was not dependent on the results of the alignments 
for (x), (0), and (0') it is possible to begin the consideration of possible A^ vacua by 
considering the alignment of (u), which may be brought to the desired (1,1,1,1) by 
action of the T generator. Then, having properly aligned (u), the VEV alignment of 
{(p) is determined uniquely. Then, the VEVs of the flavons associated with the charged 
leptons {i.e. Xi 'Pi and 0') may be brought to the desired form by application of relevant 
powers of 5, T, and U matrices. In general, these transformations will change the already 
aligned (u) by multiples of T. This in turn will force (yj) out of the alignment needed 
for GR mixing. Therefore, Nature is required to choose one of five vacua necessary for 
the Golden Model. This is a mild assumption. 

4 Phenomenological Predictions 

In this section of the paper, we study the phenomenological implications of the complex 
mass sum rule of Eqs. f|T7|l and fH3l) . Although this sum rule has been studied before j28j . 
this work was done in the era of a small reactor angle. There has been some mass sum 




(59) 
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rule studies after the advent the measurement of a large value for ^13 [29j, but this work 
did not go into great detail for the inverse mass sum rule we have generated in our 
model. Hence, we are motivated to perform our own analysis. 

We begin the analysis of the phenomenological implications of the Golden Model by 
re-expressing the complex mass sum rules of Eqs. f|T71) and (H3|) in terms of the physical 
neutrino masses, \mi\, and Majorana phases, 0:2 and 0:3 (note the sign convention change 
on Oj): 

(60) 

From the above 
(61) 



+ 

\mi\ \m2\ \m-i\ 

where the superscripts "LO" and "NLO" have been dropped on rrii 
formula, it is easy to see that Eq. f lBU]) implies that 



1 



1 



mi 



\fn2\ 



2 cos A 
|m2| |m3| 



where A = 02 — 0:3- With the realisation of this form, it becomes clear that together 
with the two experimentally measured values of the neutrino mass squared differences 

|mjp — |mjp, it is possible to numerically calculate 



Am2i and 



Am23|, where Amjj 



the individual neutrino masses for a specified value of A. 
can be found in Figure [TJ 



The results of this analysis 



0.30 



0.25- 




0.0 0.5 1.0 

COS(A) 

Figure 1: mass (eV) vs. cos(A) for the normal and inverted neutrino mass orderings 
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The plot of Figure [T] contains many interesting features. The first of which is the 
discontinuity at cos(A) ^ .5. Namely, given the experimentally measured values for 
Amfj and the neutrino sum rule from Eq. fl6Tl) no solution exists. The second interesting 
feature is that only normal ordered neutrino masses can only exist below cos(A) ^ .5, 
and inverted only above cos(A) ^ .5. This implies that varying A provides a transition 
between normal and inverted neutrino mass orderings, the justification of why there is 
no solution for neutrino masses at cos(A) ^ .5. 

The physical masses as a function of cos(A) serve as a stepping stone for numerically 
calculating the values of ^, S, and \f3\ (c.f. Eq. f lT8|) ) for a fixed value of cos(A). To 
correctly calculate these parameters to NLO, it is necessary to use the NLO masses of 
Eq. f l42|) . However, notice that these masses include an additional complex parameter 
e = |e|e*^'. We find |e| by demanding that the TB deviation parameter r (c.f. Eq. fl5^ ) 
matches the corresponding value derived from the central value of the reactor angle 
given in Ref. [10], i.e. r = .22. Then, |e| may be solved for as a function of the ^, 
S, 1/3 1 parameters and removed from the NLO masses of Eq. (142|) leaving 4 unknown 
real parameters. can be eliminated by minimising and maximising the masses with 
respect to 6^. After this, the maximised and minimised NLO masses may be expressed 
in terms of only ^, S, and \(3\. Setting these analytic forms for the maximised and 
minimised NLO neutrino masses equal to their corresponding values found in Figure [1] 
reveals the values of the parameters ^, 5, and as a function of cos(A)jl Notice that 
the parameters entering in the masses also appear in the analytic forms of the corrected 
mixing matrix elements, c.f. Eqs. fl2^ - fl5Bl) . Therefore, it is possible to numerically 
calculate the elements of the PMNS matrix and find the TB deviation parameters s and 
a of Eqs. fHO|) - (HT|) as a function of cos(A). 

4.1 NLO Corrections to s and a 

In this section, we report the results obtained from using the previously discussed 
method for determining |e|, \/3\, ^, and 6 to calculate the TB deviation parameters 
s and a to NLO. Notice that in the definitions of a and s in Eqs. (H0l) - (l4Ti) . the un- 
known Majorana phases 02 and appear, s and a are minimised and maximised with 
respect to these phases to remove their dependence on them. Then, the calculation of 
how the s and a parameters vary as a function of cos (A) is straightforward. The plots 
resulting from parameter scans ranging over values of cos(A) G [—1, 1] for the s and a 
parameters can be found in Figures [2] and [3l respectively. 

The plots in Figures [2] and [3] show some interesting features. The first feature to 
notice is that the LO GR prediction is contained in the allowed region predicted by 
the NLO Golden Model. This is due to the fact that if A[/i2e~*~ and Af/23e~* 2 are 

^Numerical instabilities in some of the solutions for the 6, and \f3\ results of the minimised NLO 
neutrino masses urged us to use the maximised results in our analysis. Although, the stable results 
show little difference between the two cases. 
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Figure 2: s vs. cos(A) for fixed r = .22: Tiie region between the red lines represents the 
possible values the s parameter can take for normal ordered neutrino masses, and the region 
between the light blue lines represents the possible values for the inverted ordering. The 
thickness of the lines takes into account the la ranges for the Am?^ used in producing the 
plot, and the grey area represents the experimentally allowed la range from Eq. ([3]). The LO 
Golden Ratio Prediction is given by the dashed line. 

purely imaginary then the corrections to s and a vanish, respectively, as can be seen 
from Eqs. (HO|) - (HT|) . The second slightly more interesting feature to notice is that the 
at cos(A) = .5 there exists no solution, a relic of the mass vs. cos(A) analysisjfl Yet, the 
most intriguing result of this analysis is that the values for s and a can actually lie within 
the la band of Eq. Clearly, the NLO corrections provide the needed adjustment to 
the Golden Model's problematic LO predictions of r = a = and s = —.09, by allowing 
them to be shifted to agree with current experimental bounds. 

Having calculated the allowed ranges which s and a can take (as a function of cos (A)), 
there still exists another plot which needs to be generated to further probe the predic- 
tions of the Golden Model, a vs. s. The plot of a vs. s will further restrict our parameter 
space of the Golden Model through the correlation of a and s. The method for gener- 
ating this plot is straightforward because all relevant parameters are known with the 
exception of the Majorana phases, a2 and as, and the phase of e, 6^. Unfortunately 
A = 02 — as does not enter in either of the forms of the a or s parameters in Eqs. 
( HU]) -f l4ip . so the method of taking the maximum and minimum values cannot be used 

^It should be noted that in the graph of a vs. cos(A), the value of a in the inverted hierarchy takes 
larger values which begin to break down the r, s, a, perturbative expansion for the PMNS matrix in a. 
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Figure 3: a vs. cos(A) for fixed r = .22: Tlie region between tlie red lines represents the 
possible values the "a" parameter can take for the normal ordered neutrino masses, and the 
region between the light blue lines represents the possible values for the inverted ordering. 
The thickness of the lines takes into account the la ranges for the Am?^- used in producing 
the plot, and the grey area represents the experimentally allowed la range from Eq. ([3]). The 
LO Golden Ratio Prediction is given by the dashed line. 

here because we want a value of a for each s. Thus, it is necessary to introduce an- 
other parameter S=a2 + as, so that the parameter A may be utilised when analysing 
how s and a vary with each other. Doing this enables a2 and to be re-expressed as 
a2 = 1/2(A + S) and as = 1/2(E — A). Since all parameters are known for each value 
of A, all that is left unknown is S and 6^. Generating random values for E,^^ G [0, 27r], 
while keeping the consistency of the definition of the Majorana phases, yields the plot 
found in Figure HI As can be seen, this plot further constrains the parameter space of 
the Golden Model, as it contains a subset of the values displayed in the plots of Figures 
[2] and [3l through the correlation of a and s. Yet even though the correlation of a and 
s further reduces the Golden Model's parameter space, the experimentally determined 
values for a and s can be still accommodated by the Golden Model, providing consistent 
phenomenological predictions which are able to match the bounds of Eq. How- 
ever, to further analyse the phenomenology of this model, we next consider the Golden 
Model's predictions for neutrinoless double beta decay. 
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Figure 4: a vs. s for fixed r = .22: Purple represents predictions from normally ordered 
neutrino masses and the light blue points from inverted neutrino masses. The purple region is 
also where they overlap. The black rectangle designates the allowed la region of Eq. Q and 
the black cross is the LO Golden Ratio prediction. 

4.2 Neutrinoless Double Beta Decay 

From the plots in Figures |2] and [3l it is seen that as the parameter A varies, a transition 
between normal ordered and inverted ordered neutrino masses occurs. Therefore, it is 
interesting to analyse the behaviour of the effective Majorana mass parameter, as 
a function of cos(A), where 

mpp = \\mi\U^^ + |m2|t/iV + |m3|f/iV""'l- (62) 

This analysis is preformed by utilising the previously discussed parameters A = 02 — «3 
and E = ^2 + as, so that Eq. (l62l) may be re-expressed in terms of S and A, instead of 
a2 and 03. Then, it is possible to invoke the work from previous sections to numerically 
generate a plot revealing how m^^ changes as cos(A) varies. The results from this 
analysis can be found in Figured Notice that (as postulated), cos(A) "disentangles" 
degenerate neutrino orderings, as can be seen from the transition at cos(A) = .5. Also, 
the results are consistent with the experimental upper bounds on m^j^ reported by EXO- 
200 [30] . Yet the most important conclusion one can draw from Figure [5] is that the 
Golden Model is testable in the near future. Due to the neutrino mass sum rule of 
Eqs. f[T7l) and fH51) . the normal neutrino mass ordering's prediction for m^^ cannot 
be arbitrarily small. Hence, the model has the capability of being tested by the next 
generation of neutrinoless double beta decay experiments. 

We conclude the analysis of the Golden Model's phenomenology by constructing, 
niisp vs. miightest- This is possible to do since the neutrino masses and m^p are known 
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Figure 5: m^^ vs. cos(A): The red region corresponds to the values that can take in 
the normal neutrino mass ordering, the light blue region corresponds in the inverted ordering, 
and the grey region is the experimental range of upper bounds that m^^g can take from EXO- 
200 [30]. 

for each value of cos(A). Thus, it is straightforward to find the maximum and minimum 
values of m^^ as well as the value of the lightest neutrino mass {i.e. mi in the normal 
ordering and in the inverted ordering) for each value of cos(A). These results can 
then be combined to produce the gold regions of the plot found in Figure [61 The plot 
in Figure |6] also contains blue and red regions corresponding to the normal and inverted 
neutrino mass orderings if there was no neutrino mass sum rule. The different shades 
of the blue and red regions correspond to the the different ranges obtained when taking 
into account the best fit and la deviations away from the best fit0 It should be clear 
from Figure [6] that the Golden Model's inverse mass sum rule severely restricts the 
allowed values that can take, leading to a lower bound of mughtest ~ 0.01 eV and 
mpp > 0.005 eV for the normal mass ordering, with the mass spectrum extending into 
the quasi-degenerate region. As previously discussed, this allows the Golden Model to 
be tested by near future neutrinoless double beta decay experiments. 



^"^It is important to note that in order to generate the model independent case of Figure we have 
used the code developed in Ref . ^Slj with the global fit values of Ref . . 
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Figure 6: mjSis vs. mughtest- The red and light red regions represent the model independent 
values that the inverted neutrino mass ordering can take based on the central value and la 
deviation of Ref. [10], respectively. The blue and light blue regions are the analogue of this 
for the normal neutrino mass ordering. The gold regions correspond to the Golden Model's 
prediction for m/3/3 in both the normal and inverted orderings. 

5 Conclusions 

The Golden Ratio mixing prediction of a zero reactor angle means that this model, 
along with other simple schemes such as TB mixing, is no longer viable. Similarly, 
there is increasing evidence that the atmospheric angle lies in the first octant, and 
so the prediction of maximal atmospheric mixing, also shared by TB mixing, is also 
disfavoured. However, while the TB prediction for the solar angle remains viable, the 
Golden Ratio solar angle prediction is also disfavoured. This provides a motivation 
for considering the Golden Ratio predictions at NLO. Unlike the case of TB mixing, 
however, where the LO solar angle prediction is viable, and so one may preserve a 
subgroup of the Klein symmetry, in the case of the Golden Ratio predictions there is a 
strong motivation to completely break the original Klein symmetry at the NLO. 

In this paper we have proposed a new model of leptons which corrects the LO 
predictions of Golden Ratio mixing via a minimal NLO Majorana mass correction which 
completely breaks the original Klein symmetry of the neutrino mass matrix. The min- 
imal nature of the NLO correction leads to a restricted and correlated range of the 
mixing angles allowing agreement within the one sigma range of recent global fits. Re- 
markably, the minimal NLO correction also preserves the LO inverse neutrino mass sum 
rule, implying a neutrino mass spectrum which extends into the quasi-degenerate region. 
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allowing the model to be accessible to the current and future neutrinoless double beta 
decay experiments. 
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Appendix A: The Group Theory of 



Ar, is a simple groupo By definition, it is the group of all even permutations of a 
set of five elements. As such, it has 5!/2 = 60 elements. In fact, it can be shown 
to be isomorphic to the Icosahedral Symmetry Group, X, the group of all rotations of 
an icosahedron which preserve the icosahedron's orientation. This isomorphism turns 
out to be insightful when considering the conjugacy classes of A^ = X, as one can 
express the conjugacy classes in terms of Schoenflies notation where they are denoted 
by and represent rotations by — . With the further definition that the number in 
front of is the number of elements contained in a conjugacy class, the conjugacy 
classes of A^ can be written as: / (the conjugacy class consisting solely of the identity), 
156*2, 2OC3, I2C5 and 12C|. Using these conjugacy classes and the theorems that 
posit that the number of elements of a group is equal to the sum of the squares of 
the dimensions of the irreducible representations and that the number of irreducible 
representations of a group is equal to the number of conjugacy classes, it is easy to 
see that A^ has a one dimensional irreducible representation (1), two three dimensional 
irreducible representations (3 and 3'), one four dimensional representation (4), and one 
five dimensional irreducible representation (5) because 

1 + 15 + 20 + 12 + 12 = 60 = 1^ + 3^ + 3^ + 42 + 5^. (A.l) 

Yet in order to do anything physically useful with these five irreducible representations, 
an explicit matrix representation for each irreducible representation of A^ must be cal- 
culated/found. To do this it is helpful to find a set of generators and rules defining their 
multiplication {i.e. a presentation) which generate ^5. We choose to work in the basis 
given in Ref. [20] in which is generated by two generators, S and T, satisfying the 
presentation rules 

S^ = T^ = (ST)^ = 1. (A.2) 

The explicit forms of S and T can be found in Ref. |20] . For completeness, we choose to 
list them again here along with another element of A5, U = ST"^ ST^ S . Notice that 
the S and U elements generate Z2 Klein subgroup of A^. Recall that in Section 

12.31 of this work, spontaneously breaking A^ to this Klein subgroup (in the neutrino 
sector) was used to generate Golden Ratio mixing. The explicit forms of S, T and U 
for each irreducible representation can be found in Table IH 

Each of the elements in Table H] will leave invariant a different fiavon VEV alignment 
(up to proportionality), i.e. G{(t>) = {(/)) where G = S,T or U and {(f)) is the VEV of 
a fiavon field preserved by the action of the group element G. These alignments are 
crucial when considering the spontaneous breaking of A^ to having fiavons 

develop VEVs invariant under the action of the S and U preserves the Klein subgroup. 
A listing of these "invariant" alignments is given in Table |5l 

^^A^ contains no nontrivial normal subgroups. 
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Table 4: The S, T, and U elements of for the 1-, 3-, 3'-, 4-, and 5-dimensional irreducible 
representations in which p = e~ and (f)g = {1 + 75) /2 is the Golden Ratio. 
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Table 5: The Invariant VEV Alignments for one and only one of the S, T, and U elements. 
The Wi and Vi represent components of the 4- and 5-dimensional flavon VEVs left unspecified 
by S and U, respectively. Notice that there exists no non-trivial VEV invariant under the 
action of T4. 

Notice that if both S and U generators are unbroken, such that the Klein symmetry 
exists at low energies, then the only nontrivial alignments that exist are for the 4- and 
5-dimensional irreducible representations. The VEV which preserves both ^'4 and U4 is 
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given by 



/ 1 \ 



1 
1 



(A.3) 



\1/ 



whereas the VEV which preserves 5*5 and is given by 



V2 
V3 



(A.4) 



\ V2 / 



Recall that the five dimensional VEV given above was crucial in constructing the LO 
Golden Model of Section 12.31 In addition to constructing these invariant VEVs the S 
and T elements/generators can also be used to calculate the Clebsch-Gordan coefficients 
associated with the particular irreducible representations of A^. This is the goal of the 
remainder of this appendix. 
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Appendix A.l: The Kronecker Products and Clebsch-Gordan 

Coefficients of 

From the matrix representations of the S and T elements/generators of given in the 
previous section, it is straightforward to calculate the Clebsch-Gordan coefficients for 
the decomposition of the product representations, which we now list for this basis in 
detail. We use to denote the elements of the first representation, bi to indicate those 
of the second representation of the product, and the subscripts "a" and "s" to indicate 
a representation which is antisymmetric or symmetric, respectively. 



3 3 Ig ^ 3^ ® 5g 



aibi + 0263 + 0362 

3a ~ I 0162 - 0261 
0361 - 0163 



f 2aibi - 0263 - 0362 \ 
— \/3ai62 — -\/3a2&i 
1/602^2 
VGasbs 
V -V^aibs - V^asbi J 



3' (g) 3' = Is ® 3; ® 5s 
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aibi + 0263 + 0362 

^2^3 - ^3^2 
0162 - 0261 
0361 - 0163 

/ 2aibi - 0263 - 0362 \ 

\^aib2 — \/3a2^i 
\/3ai63 — y/Sa^bi 
V6a2b2 



I 



3 ® 3' = 4 ® 5 



/ \pla2bi + 0362 \ 
-\pla1b2 - a-ibj, 
-1/20163 - 0262 
V \f2a-ibx + 0263 ) 



( \/3ai6i \ 
0261 - v^a362 
0162 - sf^a-ib-i 
aibs - 1/20262 
\ 0361 - V^a2b3 ) 
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l(g)4 = 3'®4®5 
— -\/2a2&4 — V^asbi 
"\/2ai63 + 0262 - 03^*4 



/ axhx - V2a3b2 \ 
—0162 — V^a2bi 
0163 + V^asbi 
\ -0164 + \/2a2bz / 



( \/6a2b4 - y/Qa^bi \ 
2\f2aibx + 20362 
—\f2axb2 + 0261 + 80363 
-\/2ai63 — 80262 ~ ^3^4 
V -2V^oi64 - 20263 / 



3^04 = 30405 

-\/20263 - V^0362 

a/2oi6i + 0264 - 0363 
"\/2oi64 - 0262 + 0361 



/ 0161 + ^20363 \ 
O162 — -\/2o364 
-0163 + "\/2o26i 

\ -O164 - V^0262 / 



/ ^60263 - ^60362 \ 
\f2a\b\ — 3a2b4 — 0363 
2V20162 + 20364 
— 2-\/2ai63 — 20261 
\ -"\/2oi64 + 0262 + 30361 / 



305 = 30 3' 0405 



3' 05 = 30 3' 0405 



3' 



v 



-2oi6i + \/3o265 + \/3o362 

v/30162 + O261 - V^0363 

V^aibs — "\/6o264 + 0361 

VSaibi + 0265 + 0362 
0163 - \pla2b2 - \PiazbA, 
aib4 - -\/2o263 - V^aab^ 

2^20162 - "\/6o26i + 0363 \ 

—\/2aib3 + 20262 — 30364 
v/20164 + 30263 - 20365 
— 2\/2oi65 — O264 + \/6o36i / 



3' 



/ V^0265 - v/30362 \ 

-O162 - V^026i - V^0363 

— 20163 — \/2a2b2 
20164 + \/2a3b^ 
0165 + \/2o264 + \/3o36i / 



\/3oi6i + 0264 + 0363 
O162 - V2a2b5 - V2a3b4 
O165 - V2a2b3 - V2a3b2 

— 2oi6i + \^a2b4 + ^/Sa-^b^ 
A/30163 + 0261 - VGa^b^ 
\/3oi64 — \f^a-^2 + ^3^1 

"\/2oi62 + 80265 - 20364 \ 

2v^oi63 — v^026i + 0365 
— 2-\/2oi64 — 0262 + -\/6o36i 
-\/2oi65 + 20263 - 30362 / 

-\/30264 — -\/30363 \ 

20162 + -\/2ci364 

— O163 — V^026i — V^0365 

0164 + \pla2b2 + ^730361 

-20165 - \/20263 / 
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4 (g) 4 = Is e 3a e 3^ e 4s e 5s 


4(g)5 = 3e3'e4e5ie52 


Is ~ aibi + 0263 + 0362 + 0461 

/ -0164 + 0263 - 0362 + 04^1 \ 
3a ~ V^a2b4 - V^a4b2 

\ V2aib3 - V2a3bi J 

/ 0164 + 0263 — 03^2 — 04^*1 \ 
3'a ~ ^20364 - V2aib3 

\ \/ Ia\b2 — w 2,a2b\ j 

I 0264 + 0363 + 0462 \ 

. a\bi + 0364 + 0463 

0162 + 0261 + 0464 
\ 0163 + 0262 + 0361 / 

/ \/3ai64 — \pia2b3 - \f^a3b2 + ^Pia^bx \ 
— \/2a264 + 2\/2a363 — \/2a/^2 
5s ~ —2\pia\b\ + \pla3b/^ + \piai,bz 
\pia\b2 + \/2a2^i — 1\plai^b\ 
\ - ^/2a 1 63 + 2 \/2a2 62 - \f^a?, h 1 


/ 2x/2ai65 - ^/2a2bi + \/2a3&3 - 2^/20462 \ 
3 '-^ — \/6ai&i + 20265 + 80364 — 0463 j 
\ 0164 — 80263 — 20362 + -\/6o46i / 

/ ^20165 + 2^/20264 - 2^/20363 - ^/20462 \ 
3' ~ j 80162 — \/6o26i — 0365 + 20464 

\ —20163 + 0262 + \/6o36i — 80465 / 

/ \/3ai6i — \/2a265 + -v/20364 — 2^/2a4b^i \ 
^ — -\/2ai62 — ■\/3a26i + 2\/2o365 + \/2aibi 
1/20163 + 2^/20262 — \/3o36i — \/2o465 
\ -2V20164 + \/2o263 - V20362 + \/3a46i / 

/ \/2oi65 - \/2o264 - \/2o363 + 1/20462 \ 
— \/2oi6i — \/3o364 — \/3a463 

51 ~ -\/3ai62 + -\/2o26i + \/3a365 

•/3a262 + •/2a36i + ^0465 
\ — \/3ai64 — \/3o263 — \/2a46i / 

/ 20165 + 40264 + 40363 + 20462 \ 
4oi6i + 2V60265 

52 ~ — \/6oi62 + 20261 — \/6a365 + 2\/6a464 

2-\/6ai63 — ■/60262 + 20361 — 'J&a^b^ 
\ 2\/6a362 +40461 / 
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5 ® 5 = Is © 3a © 3; © 4s © 4a © 5i,s © 52,s 



Is ~ CLibi + 0265 + a3b4 + 0463 + 0562 



3„~ 



a2&5 + 20364 — 20463 — 0562 
-\/3aib2 + V3a2bi + \/2asb^, — ^/2a5bs 
\/3ai65 + \/2a2b4 — ^J2ai^2 — ^Ci^bi 

20265 - 03^4 + 0463 - 20562 
\pia\bz — ^a-ib\ + \/2a,\br-^ — \f2a^b,^ 
--\/3cii64 + \pia-^-i — \/2a-^2 + \/3a46i 



/ 3-\/2cn62 + 3\/2a26i — ^f?>a■^^ + ^^f?>a^^/^ — ^f?)a^b•i \ 
3-\/2ai63 + 4-\/3a262 + 3-\/2ot36i — ^^ia^^ — -\/3a564 
3\/2oi64 — \/3o263 — a/3o362 + 'i\/2a,\b\ + 4^/30565 
\ 3-\/2ai65 — -\/3a264 + 4-\/3a363 — \fZa/^2 + 3-\/2a56i / 



/ \/2oi62 - A/20261 + V30365 - V30563 \ 
— \/2ai63 + \f2a-ib\ + \/3o465 — 1/30564 

— -\/2oi64 — -\/3ct2^3 + '\/30362 + '\/2046l 

\ -\/2ai65 — -\/3ci264 + -\/3a462 — V^osfti / 



»l,s 



/ 2oi6i + 0265 - 20364 - 20463 + 0562 \ 

0162 + O261 + \/60365 + A/6O563 

-20163 + \/6o262 - 20361 
—20164 — 20461 + -\/6a565 
\ 0165 + "\/6a264 + \f^a/^2 + a^b\ ) 



( 2oi6i - 20265 + O364 + 0463 - 20562 \ 
— 20162 — 20261 + \/60464 

0163 + 0361 + -\/6a465 + ^f^a^b/^ 

0164 + "\/6a263 + \f^azb2 + 0461 
V -20165 + v/60363 - 20561 / 



33 



Appendix B: A Basis of Contractions for AW^, 

In what follows, we present a basis for the contractions of the NLO corrections to the 
Majorana neutrino mass matrix dictated by the NLO correction to the superpotential 
AWi, = XXNN. In Eq. f lT9|) . we presented the NLO PMNS correction arising from the 
specific contractions 



AW,. 



yc 



(B.l) 



In order to justify this minimal set of contractions used, we will construct all possible 
other contractions and show that these two contractions form a basis for the set. We 
begin this by noting that there are two possible pairwise contractions which one can 
construct that contain a singlet: 



((A^A^) (AA)) and ((A^A) (A^A)) . 



(B.2) 



Then, using the fact that 3 ® 3 = Ig © 3a © 5s, there are only 6 possible contractions, 
3 for each pairwise object because the singlet is contained in the Kronecker product of 
identical irreducible representations in A^. These contractions are (labelling components 
of A as and presenting the Majorana matrix arising from each contraction): 



A=iNN),^ (AA),^ 

C={{NN),^ {XX),X^ 
D = {NX),^ {NX),^ 
E={{NX),JNX)J^^ 

F={{NX),ANX),X^ 







al + 2a2a3 





aj + 2a2a3 




- a2a3) 
Qaia2 



6ai 

-2aj + 2a2a3 



6aia2 
—2al + 20203 



6an 



/ 2 






2 

(i^d^ d^ 






.d\d2 d2d^ 


a\ ) 




^-20203 




aia2 




al 




. a\a2 — 


2 

— Ci2C^3 


al 



Aa\ + 60203 
0102 



0103 



3ai + 70203 



01O2 
3oi + 702O3 



(B.3) 
(B.4) 
(B.5) 
(B.6) 
(B.7) 
(B.8) 



al 
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From the above results, it becomes clear that A and C form a basis for the other 
contractions presented because 



B = OD + OA, 
D 



E 
F 



C A 

"6^3"' 
C _2A 

"6 ^ 

C lOA 
"6 ^ ~3~' 



Going beyond pairwise contractions and considering contractions such as 

(((^^)5.A)3A)^ , 



(B.9) 
(B.IO) 

(B.ll) 

(B.12) 

(B.13) 



gives nothing new; in this case, we simply obtain C. The full list of contractions of this 
type is as follows (there are 12, up to obvious symmetric or antisymmetric permutation): 







\ ' J-s 


= D, 


(B.14) 


(((^^)3.A)3^A; 




(((^^)3.^^)3.^),^ 


= E, 


(B.15) 


(((^^)5.A)3A^ 




\ ^ J-s 


= -F, 


(B.16) 


(((AA),^iV)3iv] 






= D, 


(B.17) 






(((^^)3.A)3^iv)^^ 


= -E, 


(B.18) 






{imh^N)^x)^y- 


= -F. 


(B.19) 



Therefore, it is sufficient to only consider contractions arising from A and C, justifying 
the original assumption of Eq. f[T^ . 



Appendix C: Breaking the GR Klein Symmetry 

In this appendix we show that the assumed correction to the superpotential (c.f. Eq. 

(DSD), 

AW^ = f (NN),^ iXX\^ + ^ {(NN),^ (AA),J^^ , (C.l) 

is sufficient to completely break the Klein symmetry associated with the unsuccessful 
LO Golden ratio prediction. We begin by first noting that so long as ai 7^ and/or 
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02 7^ — ^ non-zero value of ^13 will be obtained because a VEV proportional to 
(0, 1, — 1)^ will be preserved by 




(C.2) 







Furthermore, by construction, any representation formed from the combination of two 
triplet fields preserving U3 will preserve said representation's corresponding U. For 
example, the 5 constructed from two A fields where the are left arbitrary is 



((AA), 



/ 2aj - 20203 \ 
—2y/3aia2 

A/Gag 
V -2^/3aia3 J 



(C.3) 



In the above result, the subscript "Sg" on (AA) denotes that we are selecting the Sg from 
the tensor product 3 (g) 3 = Ig © 3a ® 5g. Since 



/ 1 \ 

1 

1 

10 

\ 1 y 



(C.4) 



it is clear to sec that only two solutions exist which leave C/5 invariant {i.e. satisfy 
t^5((AA)5^) = ((AA),;): 



Oi = and 02 = —03 
02 = 03 and Oi arbitrary 



(C.5) 



From these two solutions, the latter may be interpreted as resulting from general invari- 
ance of a VEV acted upon by C/5 and the former existing because the 5g was constructed 
from two 3- dimensional irreducible representations. In cither case, the assumed VEV of 
(A) = (0, 0, 03)"^ breaks U3 and the corresponding ((AA), ) breaks f/5 providing a correc- 
tion to the problematic vanishing reactor mixing angle. Next, we turn to the invariance 
of ((AA),^) under the S generator. 
Since 



-^3 = 



1 



( 1 

-x/2 
\^ -72 



-v^ -^2 \ 



(C.6) 
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it is straightforward to see a (A) oc (— a/20c,, 1, 1)^ is left unchanged by the action of ^3 
(c.f. Table E]). This, in turn, implies the conditions for the invariance of a triplet VEV 
under the 3- dimensional S generator: 



ai = —\pl(^gaz and 02 



as- 



(C.7) 



Then, acting 



/ -1 V6 

\V6 <Pl 



2 



V6 \ 
-20, 



(C.8) 



on ((AA)5) implies 5*5 is preserved by the ((AA)^) VEVs constructed from 
-03 V2 \ / 72(03 - 1) 

\\\ \ 



(A)oc 





(C.9) 



Notice that the first two solutions preserve 6*3 and 6*3/, respectively. The final two 
alignments preserve (when it acts on their corresponding ((AA)g) contraction). From 
these S'-preserving alignments it is easy to see (A) = (0, 0, 03)"^ will also break the 
S generator. Furthermore, it is interesting to see that, as noted in Ref. [20], it is 
impossible to preserve both 5*3 and Uj, simultaneously by any VEV given to a flavon 
field transforming under the 3-dimensional irreducible representation. However, this is 
exactly what is needed to correct the problematic predictions of LO golden ratio mixing 
given in Eq. ((E 
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